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Numerical optimization procedures have been considered for the design of airfoils in transonic 
flow based on the transonic small-disturbance (TSD) and Euler equations. A sequential approxi- 
mation optimization technique was implemented with an accurate approximation of the wave drag 
based on the Nixon’s coordinate straining approach. A modification of the Euler surface boundary 
conditions was implemented in order to efficiently compute design sensitivities without remeshing 
the grid. Two effective design procedures producing converged designs in approximately 10 global 
iterations were developed: interchanging the role of the objective function and constraint and the 
direct lift maximization with move limits which were fixed absolute values of the design variables. 
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INTRODUCTION 

Current aerodynamic design methods can be broadly categorized as inverse methods, e.g., Volpe 
and Melnik 1 , and numerical optimization methods, e.g., Vanderplaats and Hicks 2 . In general, in- 
verse methods have been widely used as design tools primarily due to their computational efficiency. 
They do have a weakness associated with the closure problem which generally requires considerable 
design experience. Also inverse methods, initially developed for potential flows, have been success- 
fully applied to rotational flows based on the Euler equations, e.g., Giles, Drela and Thompkins 3 
and the Navier-Stokes equations, e.g., Malone, Narramore and SankaT 1 . 

Numerical optimization methods have not been widely used in practical airfoil designs primarily 
due to the large amounts of computational resources needed. Nevertheless, the methods will con- 
tinue to be developed since they have many advantages such as automated design capability, ability 
to handle multi-point design and varieties of constraints along with a capability of inclusion into 
multi-disciplinary design of complete vehicles. A major reason for the large computational effort of 
numerical optimization methods is the very large number of transonic analyses needed to develop 
converged designs. Some improvements to the efficiency of numerical optimization methods have 
been obtained through the implementation of the shape functions, by Vanderplaats and Hicks 2 and 
Aidala, Davis, and Mason 5 , and through the use of efficient optimization procedures, Vanderplaats 6 , 
Joh. Grossman and Haftka 7 and Joh 8 . 

The motivation for the present work stems from plans to incorporate transonic airfoil designs 
within an integrated aerodynamic/structural design of an aircraft wing, e.g., Grossman et al. 9 . 
Thus, our objective is to develop efficient numerical optimization procedures for the design of two- 
dimensional airfoils at transonic speeds, using as few complete transonic analyses as possible. 

A preliminary study for this effort is reported in Ref. 7, where some special treatments were 
developed for design optimization based upon the transonic small-disturbance (TSD) equations. In 
this paper, we will amplify and improve these ideas and examine in detail the applicability of the 
methods to the more accurate Euler equation analysis. 

DESIGN FORMULATION 


Design Problem 

The design problem considered can be stated as: 

maximize Q(X) 

such that C d {X) < C dl , (1) 

A(X) > A mtn> 

where X is the vector of design parameters X = (Xi, X 2 ,...,X N ) T specifying the airfoil geometry, 
C d is the drag coefficient due to wave drag. C dl is the prescribed upper limit on wave drag, A is 
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ANALYTICAL FORMULATION 


TSD Analysis 

The first approach taken is based on a TSD formulation. This approach is based on an asymp- 
totic expansion of the perturbation potential in the limits of airfoil thickness (f/c) — <■ 0 and free- 
stream Mach number M ’ TO — > 1, e.g.. Ref. 12. The governing equation may be written as 


(1 - Ml)<p x - 


"t" <Pyy 0 , 


with the surface boundary condition on the airfoil surface, y =Y(x) expanded to be 


(3) 


d<t>, n . dY 

ai {x ' 0) =&’ 


0 < x/c < 1 . 


(4) 


It may be noted that the governing equation retains the important non-llnear effects of transonic 
flows but does not include the effects of entropy change across shock waves. Also, since the airfoil 
surface boundary condition is applied along the axis and not on the actual surface of the airfoil, the 
calculation can be performed on a simple Cartesian grid. This is useful in the design problem, since 
we can efficiently utilize a fixed grid even with changing airfoil geometries. The specific analysis 
code used for our TSD calculations is TSFOIL described in Ref. 13. 


Euler Analysis 


In order to evaluate whether the design procedures developed using the approximate TSD anal- 
ysis are valid for more accurate flow field methods, we have investigated the airfoil design with the 
Euler equations. We utilize the complete set of governing equations for an inviscid, non-heat con- 
ducting, adiabatic flow with negligible body forces. The equations may be written in conservation-law 
form in Cartesian coordinates as 


dQ dF dC 

dt dx dy ' 


(5) 
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( 6 ) 


with velocity components u, v, density p. total energy per unit mass eo = e + (u 2 + v 2 )/2, with e being 
the internal energy per unit mass and pressure p, which for a perfect gas may be expressed as 
p = ( 7 - l)[peo - p(u 2 + u 2 )/2] . The surface boundary conditions for the Euler equations, representing 
no flow through the solid surface may be expressed as 


v(x, Y) = Y’(x)u(x, Y) , 0 < x/c < 1 , (7) 

where a prime denotes differentiation with respect to x. Thus we see that the boundaiy conditions 
must be applied on the actual airfoil geometry, requiring a new mesh to be generated at each 
stage of the design process. We partially alleviate the computational burden of re-creating the grid 
for each geometry by assuming that design changes proceed slowly, and for a specified number 
of cycles consider the grid to be fixed to a baseline airfoil geometry. Then the surface boundary 
condition must be altered to allow a small amount of mass transpiration through the surface to 
approximately account for the changing geometiy. This procedure fits in well with the sequential 
approximate optimization algorithm used in the design process, which imposes move limits on the 
design. For each optimization cycle a baseline geometry will define the grid and the grid will be kept 
fixed throughout the approximate optimization cycle. This greatly reduces the computational effort 
for the Euler designs. 

Let the subscript 6 refer to the body surface of the new airfoil and the subscript 6 0 refer to the 
body surface of the baseline airfoil. The exact surface boundary condition on the new airfoil surface 
is the vanishing of the normal velocity V b ■ hb = 0 or, as in Eq. (7) 


t* = Y'u b . 


( 8 ) 
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the airfoil cross-sectional area, non-dimensionalized by c 2 , with c the airfoil chord and A min is the 
minimum required area. The design is performed at a free-stream Mach number of Moo = 0.75 and 
zero angle of attack, a = 0. 

This type of problem has been solved first by Vanderplaats and Hicks 2 with a full potential 
code, requiring 70 exact analyses. In Ref. 6. the same problem required 44 exact analyses with a 
sequential optimization technique and data base approach where all the previous design information 
is stored and reused for constructing higher-order approximations. 

Here, we will examine the effects of utilizing two different codes, one a more approximate tran- 
sonic small-disturbance (TSD) analysis and the second, a more exact inviscid Euler analysis. How- 
ever. due to the different approximations in the two analyses, particularly the neglect of entropy 
Jumps across the shock waves, the wave drag values are found to be different, with the TSD result 
at a lower level. In order to develop somewhat similar designs between the Euler and TSD methods 
it was found necessary to utilize a larger value of C dl in Eq. (1) for the Euler designs. 

Shape functions and Design Variables 

We have chosen to design the airfoil using shape functions following the successful implemen- 
tation of Vanderplaats and Hicks 2 ; 

N + 2 

= < 2 > 

1=1 

where Y = y/c with y being the airfoil ordinate and c the airfoil chord length. The specified 
shape functions Y t are functions of the non-dimensional abscissa x/c and the parameters X t are 
the design variables. For the shape functions here, we selected four existing airfoils (N = 4), 
namely, NACA 2412, NACA &4 x -4\2, NACA 65 2 -415 and NACA 64 2 A215. There are two additional 
shape functions for imposing the boundary conditions at the trailing edge of the airfoil. These are 
Yk+i = -i -x/c on the upper surface and zero on the lower surface, and Y n+ i = -x/c on the lower sur- 
face and zero on the upper surface. Usually with TSD analyses an open trailing edge is considered; 
here we specify this thickness to be ,0025c. For Euler analyses a closed trailing edge is utilized. 
This fixes the values of the coefficients Xs+i and Xn +2 in terms of X x , . . . , X^. 

Approximate Optimization 

When a design optimization is coupled with expensive numerical analysis code, most of the cost 
of the optimization is associated with the exact analyses and sensitivity calculations. Even with 
the most efficient transonic flow analysis code, the cost of the design process may be prohibitive if 
the analysis code and an optimization algorithm are linked together directly, so that full analyses 
are made for all the function evaluations during the design process. Instead we utilize a sequential 
approximate optimization algorithm 10 . This approach replaces the original objective function and 
constraints with approximations based upon nominal values and derivatives at an initial point. Ad- 
ditionally. move limits are used to prevent the design from moving outside the bound of validity of 
the approximations. Each approximate optimization problem is solved until an optimum is found, 
and then a new approximation is constructed there, and the design optimization process is repeated 
until convergence is achieved. An approximate optimization is typically referred to as an optimiza- 
tion cycle, and this is also the terminology used here. A key part of implementing a sequential 
approximation algorithm involves the approximation of the objective function and constraints. We 
have found that these approximations play a crucial role in the design process. The procedures that 
we have developed for approximating the lift and drag appear in detail later in this paper. 

The specific optimizer used for our study is the general purpose optimization program NEW- 
SUMT -A 11 which is based on a quadratic extended interior penalty function and Newton method for 
unconstrained minimization. The program provides the user with several approximation-switching 
and move limit strategies. 
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Now we utilize the grid system which was generated for the baseline airfoil to analyze the flow over 
the new airfoil. Thus it Is necessary to evaluate the new normal and tangential velocities on the 
baseline airfoil surface Y 0 . (The normal velocity on the baseline airfoil V 6o • will not be zero). This 
can be achieved by expanding the Cartesian velocity components at the baseline airfoil surface in 
terms of the coordinates of the new airfoil surface as 



. =t* + 0(AF), 

(9o) 


■ = u b + 0{AY ) , 

(9.6) 


where AT = Y - Y 0 . Thus using the boundary condition (8) we have v bo = u h Y' + O(AY) and 
W6„ - u b + O(AY). We can write these in terms of normal and tangential velocity components to the 
baseline geometry as 


VL = 


v ** = 


v bo -u bo YJ 

y/ 

Ubo ~ VboYp 

\A+^ 


( 10 ) 

(11) 


Using Eqs. (8)— (1 1), we can write an expression for V nio in terms of which takes into account 
the vanishing of the normal velocity at the new airfoil surface, V nb = 0, as 


VL - 


0 


1 + Y’Y' 


7 [Y' — Yq) , 


(12) 


where terms of 0(AY) have been neglected. 

The boundary conditions for the Euler calculation may be evaluated by extrapolating from the 
field points to the surface y = Y 0 values of p^. p bo and and v^. The tangential velocity V t is 
computed from Eq. (11) and the normal velocity V nto from Eq. (12). Values of the total energy are 
computed from 


(peo)fco = 


Pbo Pbo( V n ka + 
7-1 2 


(13) 


The specific analysis code used for our Euler calculations is FLOMG which is based on Jameson's 
time-stepping 14 and multigrid algorithms 15 . Although the code was developed by Swanson and 
Turkel 16 for solving the Navier-Stokes equations, we will utilize it only in the inviscid. Euler solver 
mode. 


LIFT AND DRAG APPROXIMATIONS 

We first considered simple linear approximations for the lift and wave drag, as: 


/ = /° + £ (*, - *,°) , 


(14) 


where / is either the lift or the wave drag and /° is evaluated with the initial design parameters 
The sensitivity derivatives, df/dX x are evaluated using one-sided finite-difference ex- 
pressions. The effectiveness of the lift approximation is indicated in Tables la and lb for both TSD 
and Euler analysis methods. We see that when the design variables are changed by as much as 2% . 
the linear approximation of the lift coefficient remains within a 2% accuracy, compared to the exact 
analysis, for the TSD solutions and to within 0.5% for the Euler solutions. However, this situation 
is not repeated for the drag approximation. As seen in Tables la and lb the linear approximation 
for the wave drag does not correlate closely with the exact results, with errors of 25% for the TSD 
solutions and 15% for the Euler solutions, when the design variables change by 2%. This result 


449 


Third Iniemauonal Conference on Inverse Design Concepts and Optimization in Engineering Sciences 
nCIDES-nn. Editor: G.S. Dulikravich. Washineton D.C.. October 23-25. 1991. . 

was not unexpected since transonic flows are very nonlinear and shock-wave movement must play 
an important role in the wave-drag approximation. 

In order to better approximate the wave drag, we considered the method of strained coordinates 
for perturbations of transonic flows with shock waves, introduced by Nixon 17 . This method has 
been applied for airfoil approximations by Stahara 18 . In this method, the perturbations are made in 
a strained, coordinated system where the shock remains fixed. Coordinate straining uses the axial 
position of the shock wave on the airfoil, x s and the height of the tip of the shock wave above the 
airfoil. y t ■ When x s changes by Ax 3 and y ( changes by Aj/ t , the coordinates (x,y) of each point in 
the flow field are changed by Ax and Ay given by 

Ax = s(x)Ax 3 , Ay = — Ay t ■ (15a, 6) 

Ut 


For the TSD analysis, following Ref. 17, the straining function s(x) is taken as 


s(x) = 


x(c — x) 
x 3 (c- x 3 ) ’ 


(16) 


which is valid for 0 < x < c and is equal to zero otherwise. 

For the TSD analysis, the strained coordinates are used first in the process of calculating finite- 
difference sensitivity derivatives according to the following steps: 

1 . The i-th design variable is perturbed by AX, and the corresponding Ax 3 and A y t are calcu- 
lated from a new solution of the flow field and used to approximate the shock sensitivities 

dx^Ax, dyt Ayt_ ( 17o M 

dX, ~ AX, ' dXi AX, ' 

2. The flow-field sensitivities are approximated, using Ax and Ay from Eqs. (15a,b). in terms 
of the axial velocity for the nominal design u° and the axial velocity for the perturbed flow 
field u as 

— rr(x,y) = -^rr[u(x + Ax,y + Ay) - u°(x,y)] . (18) 

us\ j AAj 

The strained coordinates are used again in approximating the flow at a new design point X as 
follows: First the new shock location and shock-tip position are calculated by a linear approximation 




(19a, 6) 


t=l ---• »=i 

Values of Ax and Ay are calculated from Eqs. ( 15a,b) and then the axial velocity it is estimated from 

4 


u(x + Ax, y + Ay) = u°(x, y) + ^ ~ ( X i ~ x ?) ■ 


i=i 


dX, 


( 20 ) 


The wave drag coefficient is determined from a contour integral of the jump in pressure across 
the shock, which may be written as 13 : 

= _^/3 M -3/4(l±l) f [u] 3 dy , (21) 

1 J shock 

where is the free-stream Mach number, 6 is the nominal airfoil thickness, 7 is the ratio of specific 
heats and [u] is the jump in u across the shock. 

For the Euler analysis a different implementation of the method of the strained coordinates is 
utilized. First, since the wave drag is calculated by integrating surface pressures, only the solution at 
the body surface needs to be approximated. The coordinate straining in y direction is not necessary. 
Furthermore, accurate solutions for the Euler wave drag were obtained' by approximating the surface 
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pressures directly, rather than computing from the flow field. Thus the surface pressure sensitivity 
and approximation, respectively were calculated as 


^-[x, y o(*)] = 2^7 (p[ x + Ax > *o(* + Ax )l “ P°l x . r o(x)]) , 


p[x + Ax, y 0 (x + Ax)] = p°[x, Vo(^)I + ^2 


r* dp[x, r 0 (x)] 


»=i 


dX, 


(Xi-Xi), 


( 22 ) 

(23) 


where Ax is the amount of coordinate straining and the superscript 0 corresponds to the value for 
the nominal design. The pressure on the perturbed geometry can be approximated by using Taylor 
series expansion about this pressure on the baseline geometry and then the wave drag coefficient 
is determined from the integration of the approximate surface pressure force in x-direction. 

The non-uniqueness of the straining function has been discussed in Refs. 17 and 18. Here, 
several tests were made to select the best straining function in terms of approximation accuracy for 
the Euler analysis. The linear piecewise continuous straining appeared to be most accurate and 
well-behaved. This has been also pointed out with the full potential method in Ref. 18. The linear 
piecewise continuous function is given by 


s(x) = 


x/x 3 ; 0 < x < x , , 

(c-z)/(c- X,) : X, < x < c, 


(24) 


was used, where x 3 is the location of shock wave. 

The results of Nixon’s coordinate-straining approximation on the drag coefficient are tabulated 
in Tables la and lb. It is seen to significantly improve the wave drag approximation, with the errors 
reducing to less than 3% for the TSD solutions and 7% for the Euler solutions, when the design 
variables change by 2%. The effect of this approximation on the airfoil pressure distributions is 
shown in Figs, la and lb. 


DESIGN RESULTS 

In this section we consider several optimization strategies for the transonic airfoil design problem 
of maximizing lift with constraints on wave-drag and airfoil cross-sectional area as given by Eq. (1). 
The minimum non-dimensional area is taken to be A min = 0.075. The wave drag constraint is taken 
to be Cd, = 0.004 for the TSD designs and Cd, = 0.010 for the Euler designs. 

Designs based on the TSD analysis 

Strategy A: Approximate Optimization with Tight Move Limits - The first optimization strategy that we 
employed consisted of imposing tight move limits in the approximate optimization procedure. The 
results of applying this strategy with two different initial designs are tabulated in Tables 2a and 2b. 
In the first case, we imposed 5% move limits in order to keep the error in the drag approximation to 
within 10%. The solution ceased to improve after 27 iterations and the move limits were tightened 
to 2.5%. At 60 iterations the move limits were further reduced to 1.25% and the solution was con- 
sidered to be converged. However, in the second case, corresponding to different initial conditions, 
employing a similar strategy resulted in a completely different design, as can be seen in Table 2b. 

In order to examine whether these solutions were local maxima we considered the following: 
denoting the first converged solution as X 1 and the second as X 2 , we defined an intermediate 
design state 

X = X 2 + {(X 1 -X 2 ), (25) 

where the parameter f may be considered the proportional distance of the intermediate design 
between X 2 and X 1 . We computed both approximate and exact value for the lift and wave-drag 
coefficients for intermediate designs with 0 < f < 1. The results are plotted in Fig. 2. From the 
upper chart in Fig. 2 we see that the lift coefficient is well-behaved between design state 2 and 
state 1. exhibiting no local maxima or minima. In the lower chart, we see that the culprit is the 
drag-coefficient constraint, which exhibits a very wavy behavior about Cd = 0.0040. Thus we can see 
that if we are a design state 2, the optimizer would prevent you from moving toward state 1, since 
that would be a direction of increasing drag. 
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Strategy B: Drag Minimization followed by Lift Maximization - The strategy that we found to be useful 
to avoid the problems associated with the drag constraint consisted of interchanging the role of the 
drag constraint with the lift objective function: We had noticed that design problems of minimizing 
drag with a constraint on the lift were well behaved. In order to solve the design problem formulated 
in this paper, we adopted the following strategy: first the wave-drag coefficient was minimized with 
a constraint on the lift coefficient of C, > 0.5 and cross -secUonal area ratio A > 0.075; then, when C d 
was below 0.004 we reverted back to the original design formulation of maximizing C\ with constraints 
of C d < 0.004 and A > 0.075. 

The results of this strategy were very good. We found that during the drag minimization phase, 
very large move limits, as large as 500% could be used without any adverse effects. During the 
lift maximization phase, move limits of 20% were imposed. This strategy seemed to be robust, and 
the solutions did not depend upon the initial data. The design history for the first case, starting 
from X° = (1 0, 0.0, 0.0, 0.0) T is tabulated in Table 3a. The design results for four different initial 
designs are summarized in Table 3b. In all cases convergence to nearly the same design result 
was obtained in 8-13 design cycles. The pressure coefficient and corresponding airfoil shape of 
the initial and final designs tabulated in Table 3a, are shown in Fig. 3. The design appears to be 
physically reasonable, with a weakened shock wave and lift increased through aft camber. 

Strategy C: Approximate Optimization with Absolute Move Limits - After we obtained successful design 
results using the strategy B. we carefully investigated the TSD solutions to determine the cause of 
the noisy drag calculation. We found that the spline interpolation routine in the TSFOIL program 
generated an irregular airfoil leading-edge geometry. Although this should not be important in the 
TSD solutions, which lose their validity at the leading edge, it clearly affected wave drag calculations 
and generated noise. We replaced the original routine with a more effective interpolation based on 
the approximate arc-length of the airfoil with a periodic boundary condition and consequently was 

able to generate fairly smooth and round noses. 

Next we attempted to directly maximize the lift with tight move limits using the new geometry 
interpolations. The design process behaved much better due to the considerably reduced noise in 
the wave drag, even though the noise was not removed completely. We felt that we should increase 
the move limits in order to get faster convergence. We also found that we could produce reasonably 
efficient designs using move limits which were fixed absolute values rather than percentages of the 
design variables. By several tests, we found that initially 0.5 could be used without any adverse 
effect and then it was reduced by half when the design did not make any improvement. The design 
results using this procedure are tabulated in Tables 4a and 4b. We experienced some convergence 
difficulties with TSFOIL using this approach. We will examine this strategy in more detail with the 
Euler analysis method. 

For the TSD designs, each exact airfoil analysis using the program TSFOIL required 10-15 
CPU seconds on the IBM 3090 at V. P. I. & S. U.. with N+l(5) analyses needed per design cycle. 
The approximate optimization using the program NEWSUMT-A required 10-12 CPU seconds on the 
same computer. 

Designs based on the Euler analysis 

On the basis of the TSD design experience, we applied the two successful design strategies. B 
and C to the same design problem with the Euler analysis. Recall that the wave-drag constraint 
value was changed to 0.01 for the Euler design due to the differences in wave drag prediction between 
the TSD and Euler methods. The original value of 0.004 was found to be too stringent for the design 
problem with the Euler analysis. 

Strategy B • Drag Minimization followed by Lift Maximization - Table 5a represents the complete design 
history for the first case, starting from X° = (1.0, 0.0, 0.0, 0.0) T . In the lift maximization phase, 
initially 50% move limits were utilized, which yielded a large improvement in the lift coefficient to 
a value of C t = 0.7136 for the first 8 iterations. After that we reduced move limits by half twice and 
then finally we imposed very tight move limits of 2% which after 25 iterations resulted in Q = 0.7144. 
For the purposes of this study, we consider the design achieved after 8 iterations, corresponding to 
the 50% move limits to be acceptable as a final design. The equivalent value of C, in the table is 
the relevant estimated lift coefficient when all of the violated constraints are brought to be critical. 
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Lagrange multipliers are used to estimate these equivalent lift coefficients since they represent the 
amount of change In the objective function due to the unit change in a constraint at an optimum. 
This provides us information whether the design makes true improvement or not. We found that 
this strategy was as robust and efficient for the design based on the Euler method as it was for the 
TSD method. As shown in Table 5b. four different initial conditions yielded nearly the same final 
design in 10-12 design cycles. The pressure distributions and corresponding airfoil shapes of the 
initial and final designs are shown in Fig. 4. 

Strategy C: Lift Maximization with Absolute Move Limits - The complete design history for the first 
case with strategy C is tabulated in Table 6a. Here also, we consider the design achieved after 
9 iterations, corresponding to the 0.5 move limits to be acceptable as a final design. Table 6b 
represents the summary of the design results for four different cases. All cases did converge to 
approximately the same design result in S-10 global iterations. 

For the Euler designs, each exact airfoil analysis using the program FLOMG required approx- 
imately 40 CPU seconds on the CRAY 2S at NASA Langley. The approximate optimization using 
the program NEWSUMT-A required 40-60 CPU seconds on the IBM 3090 at V. P. I. & S. U. The 
additional computer time associated with the approximate optimization of the Euler design is partly 
related to the more complicated wave drag calculation compared to that used for the TSD design. 

Error Magnification during Optimization 

Table 7 compares the lift/wave drag ratios predicted by TSD and Euler methods for the four 
airfoils used in the shape definition and the optimum TSD airfoil. It is seen that the agreement 
between TSD and Euler is much poorer for the optimized airfoil. This indicates that there may be 
a risk associated with optimization based on an approximate method. The optimization procedure 
may “improve" the design by exploiting the weaknesses of the approximation. 

CONCLUSIONS 

We have considered numerical optimization procedures for the design of transonic airfoils based 
on the transonic small-disturbance (TSD) and Euler equations. A sequential approximate optimiza- 
tion procedure was implemented with accurate approximation of the wave drag based on the Nixon’s 
coordinate straining technique. A modification of the surface boundary conditions was utilized in 
order to efficiently compute sensitivity derivatives without remeshing the grid with the Euler anal- 
ysis. 

The airfoil design problem which we considered consisted of maximizing the lift with constraints 
on the wave drag and area. We found that when the computed drag did not vary smoothly with the 
design parameters, the optimization process produced local extrema. A procedure interchanging 
the role of the objective function and constraint, initially minimizing drag with a constraint on the 
lift was found to be effective in producing converged designs. This procedure was also proven to be 
robust and efficient for cases where the drag varied smoothly, such as with the Euler solutions. The 
direct lift maximization with move limits which were fixed absolute values of the design variables, 
was also found to be a reliable and efficient procedure for designs based upon the Euler equations! 
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change* 

% 

Lift Coeff. 

Drag Coeff. 

Linear 

Exact 


Coord- St. 


0.0 






0.5 

0.5472 

0.5471 

0.00584 

0.00579 

0.00585 


0.5210 

0.5216 

0.00439 

0.00447 

0.00436 

1.0 

0.5602 

0.5614 

0.00656 

0.00654 

0.00657 

-1.0 

0.5080 

0.5098 

0.00367 

0.00396^ 

0.00400 

1.5 



0.00728 

0.00736 


-1.5 

wtm 


0.00295 

0.00346 

0.00342 

2.0 

0.5864 

0.5918 

0.00800 


0.00845 

\mm 

0 4818 

0.4881 

0.00223 

0.00302 

0.00295 


* Design variables increased by specified percentage from 
(0.5, 0.5, -0.5, 0.5) r 

Table la. Lift and Drag Approximations - TSD. 


change* 
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Lift Coeff. 1 

Drag Coeff. 

Linear 

Exact 

Linear 

Coord-St. 

Exact 1 

0.0 






0.5 

0.4948 

0.4942 

0.00866 

0.00870 

0.00854 


0.4809 

0.4813 

0.00708 




0.5016 

0.5006 

0.00946 

0.00957 

0.00923 

-1.0 


0.4749 

0.00628 

0.00652 

0.00664 

1.5 

0.5086 

0.5070 

0.01029 

0.01049 


-1.5 

0.4667 

0.4685 

0.00549 

0.00590 

0.00607 

2.0 

0.5153 

0.5135 

0.01104 

0.01149 


-2.0 

0.4598 

04621 

0.00470 

0.00533 

0.00553 1 


* Design variables increased by specified percentage from 

(0.5, 0.5, -0.5, 0.5) r 

Table lb. Lift and Drag Approximations - Euler 
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Table 2b. TSD Design Strategy A : Approximate Optimization 
with Tight Move Limits - initial condition 2. 
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cycle 

Design Parameters 
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mm 


mm 

mm 

0 

0.8 

-0.4 

0.7 
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0.0739 

1 

mm 
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0.0038 

0.0729 

- 

— 

— 

— 

— 


ijfM 

n ■ 

9 

1.093 

-0.483 

0.504 

-0.273 




- 

— 

— 

— 

— 

— 

— 

■ 

11 

1.114 

-0.487 
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-0.275 

0.4049 

0.0040 


- 

— 

— 

— 

— 

— 

— 

HUM 

15 

1.122 

-0.487 

0.487 

-0.277 

0.4055 

0.0040 
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cycle 

Design Parameters 

Ci 
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mm 

X, 

mm 

mm 

0 



-0.5 

0.5 

0.5296 


0.0805 

1 
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0.4954 


0.0779 

- 

ms 

im 

gpf 

— 

— 


— 

27 




0.481 

0.5385 


0.0750 

- 


bh 

JPP 

- 

— 


— 

60 







0.0750 

- 


gpl 


Sim 

■ 

nil 

- 

71 



-0.805 

0.591 

0.5592 

0.0040 

0.0750 


Table 2a. TSD Design Strategy A : Approximate Optimization 
with Tight Move Limits - initial condition 1. 
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Design 

cycle 

Design Parameters 

Ci 

H 

■ 

mm 

■9 

mm 

mm 

0 


0.0 

0.0 

0.0 

0.5656 

0.0103 

0.0822 

1 




-0.096 

0.5195 

0.0054 

0.0750 

2 

0.773 

0.280 

-0.244 

0.137 

0.5041 


0.0750 

3 




0.295 

0.4999 

0.0032 

0.0750 

4 

0.504 

0.516 

-0.476 

0.404 

0.4989 

iiimwa 


5 




0.542 

0.5005 

0.0025 

0.0750 

6 




0.945 

MKtflU 

0.0017 

0.0750 

7 

0.080 

1.308 

-1.614 

1.261 

0.5124 

0.0011 

0.0750 

8 

0.001 

1.475 

-1.938 

1.516 

0.5049 

hhhhm 

0.0750 

*9 

0.001 

1.770 

-1.863 

1.215 

0.6696 

0.0034 

0.0750 

10 




1.209 

0.7078 

0.0041 

0.0750 


* begin lift maximization. 

Table 3a. TSD Design Strategy B : Drag Minimization followed 
by Lift Maximization - initial condition 1. 


case 

Design Parameters 

Ci 

Cd 

A 

design 

cycles 


X 2 

*3 

x 4 

1 

1.0 

0.0 

0,0 

0.0 

0.5656 

0.0103 

0.0822 

10 

0.001 

1.876 

-1.940 

1.209 

0.7078 

0.0041 

0.0750 

2 

0.0 

1.0 

0.0 

0.0 

1.0676 

0.0459 

0.0771 
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-0.025 

1.878 

-1.889 

1.179 

0.7149 

0.0042 

0.0750 
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0.0 

0.8 

0.0 

0.6491 

0.0297 

0.0786 

13 

0.043 

1.904 

-2.068 

1.279 

0.7034: 

0.0041 

0.0750 

4 

0.0 

0.0 

0.0 

1.0 

0.5534 

0.0242 

0.0996 

10 

-0.006 

1.923 

-2.033 

1.270 

0.7091 

0.0040 

0.0750 


Table 3b. TSD Design Strategy B : Drag Minimization followed 
by Lift Maximization - summary of designs with var- 
ious initial conditions. 
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Design Parameters 

mm 


mm 

cycle 

mm 

mm 
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mm 

8 B 8 B 
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0 

1.0 

0.0 

0.0 
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0.0061 
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2 
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3 
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-1.249 
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0.929 
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7 



-1.504 

1.024 
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8 

-0.224 

1.916 

-1.752 
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Table 4a. TSD Design Strategy C : Approximate Optimization 
with Absolute Move Limits - initial condition 1. 


case 

Design Parameters 

c t 
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mm 

mm 

mm 
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0.0 

0.5656 



8 
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1.916 

-1.752 
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0.7218 
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1.0 



1.0676 

0.0459 

0.0771 

8 





0.7360 

0.0040 

0.0750 

3 



0.8 


0.6491 

0.0297 

0.0786 

8 

-0.539 

2.340 

-2.069 
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0.7922 

0.0041 

0.0753 

4 

0.0 




0.5534 

0.0242 

0.0996 

5 

-0.434 


-1.749 

1.224 



0.0750 


Table 4b. TSD Design Strategy C : Approximate Optimiza 
tion with Absolute Move limits - summary of designs 
with various initial conditions. 
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• begin lift maximization 
& itarting point for reduced move limits. 

** C/, f . : Equivalent Value of C\. 

Table 5a. Euler Design Strategy B : Drag Minimization followed by Lift 
Maximization - initial condition 1. 


Move 

Limit# 

Design 

cycle 

Design Parameters 

m 

c. 

■ 

H 

X , 

El! 

x* 

x > 


0 

1.0 

0.0 

0.0 

■Hf 



nimul 

1 

1 

0.652 

0.500 

-0.220 

gT%H 

.56529 

BHB 

EEJ 

2 

0.239 

1.000 

-0 672 


I.'.ilttl 


eh 

El 

g1»lr*i 

1.267 

-0.923 

Wk^zi 

F.UMI 

irmH 

eh 



1.470 

-1.098 

0.772 

64813 

BESS 


5 

-0.227 

1.635 

-1.227 


Em 

lonWjfi 


6 

-0.341 

1.792 

-1.380 

ED 

w\m 

E 33 


7 

-0.434 

1 956 

BE 

nBI 

Ers 

ebe 


8 

-0.514 

2.082 




ee^i 


*9 

-0 605 

2.159 





[|^E1 


mm 

EEM 

gg*l 







0.25 

El 

fcliU 


Blz>i 

gi&?53 






EE51 


BE 

giro 


yrmni 


■lrftM:n 



2.340 

-1.907 

wim 


ueib 


WrOlr^l 

HO 

bsh 


BEI 

gym 

IKT3EI 



SS2I 

0.125 

mm 

BSSl 

gsm 


mW'Mit 

BEH 

EuHIS 

E^3 

■IrOW^ 

ma 





ma 




0.02 




gwm 


Eu3 

yrTvTn 

eh 

E33 


- 

- 

- 

- 

- 

- 

- 




2.415 

-2.005 

1.521 


Iirarari 


BE 


- 

- 

- 

- 

- 

- 

- 

- 

35 

-0.804 

2.435 

-2.028 

ggm 





36 

-0.807 

2,438 

-2.032 

1.543 

71343 

.010000 

0.0750 

.71339 


A starting point for reduced move limits. 

** Cit f > : Equivalent Value of C<. 

Table 6a. Euler Design Strategy C : Approximate Optimization with Ab- 
solute Move Limit# - initial condition 1. 
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mm 
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1.0 
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0.0 
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11 
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3 

0.0 
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12 
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4 
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0.0 

1.0 

0.3558 

0.0165 

0.0984 

10 




1.322 

0.7090 
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Tfcble 5b. Euler Design Strategy B : Drag Minimization fol- 
lowed by Lift maximiation - summary of designs 
with various initial conditions. 
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Table 6b. Euler Design Strategy C : Approximate Optimiza- 
tion with Absolute Move limits - summary of designs 
with various initial conditions. 


Airfoils 


TSD Anavsis 

Euler Analysis 

Airfoil 1 

54.9 

41.4 

Airfoil 2 

23.3 

26.2 

Airfoil 3 

21.9 

23.9 

Airfoil 4 

22.9 

21.6 

Optimized, TSD 

171.0 

83.4 


Table 7. Error Magnification during Optimization. 



Bose solution, X = (0.50. 0 50.-0.50. 0.50) 

o-o Approximation, X = (0.51, 0.51,-0.49, 0.51) 

Exact Anolysis, X = (0.51. 0.51.-0,49. 0.51) 

Figure la. Approximation of Pressure Distribution Based On 
the Method of Strained Coordinates - TSD. 













































































































Third International Conference on Inverse Design Concepts and Optimization in Engineering Sciences 






Boss solution, X = (0.50, 0.50. 0 50, 0.50) 

0-0 Approximation, X = (0.51. 0.51,-0 49. 0.51) 
Exacl Analysis. X = (0.51. 0.51.-0.49. 0.5l) 


Figure lb. Approximation of Pressure Distribution Baaed On 
the Method of Strained Coordinates - Euler. 
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Behavior for Two TSD Local Optima 
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Figure 3. TSD Design Results 


Figure 4. Euler Design Results 









